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The zero temperature effective equation of motion is derived for a scalar field interacting with 
other fields. For a broad range of cases, involving interaction with as few as one or two fields, 
dissipative regimes are found for the scalar field system. The zero temperature limit constitutes a 
baseline effect that will be prevalent in any general statistical state. Thus, the results found here 
provide strong evidence that dissipation is the norm not exception for an interacting scalar field 
■ system. For application to inflationary cosmology, this provides convincing evidence that warm 

inflation could be a natural dynamics once proper treatment of interactions is done. The results 
found here also may have applicability to entropy production during the chiral phase transition in 
heavy ion collision. 
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I. INTRODUCTION 

Statistical mechanics generally expects an interacting system to equally distribute its energy amongst all the con- 
stituent degrees of freedom. Efforts to understand this most basic statement from quantum field theory have focused 
on simple dissipative systems. The most common example amongst these is the case of a quantum scalar field which is 
interacting with other fields and which has an amplitude that is initially displaced from equilibrium. There has been 
considerable interest in this kind of problem due to its many applications, ranging from the dynamics of condensates 
in condensed matter physics to the evolution of a inflaton field in inflationary cosmology (for a recent review see pi). 
Several methods have been applied to this problem involving analytic and numerical analysis as well as lattice 

simulations |l5|Jl6|] . 

Among the variety of dissipative dynamics that can be studied for such scalar field systems, the most tractable, 
given the present state of understanding in quantum statistical mechanics, is the adiabatic regime, where the motion 
of the scalar field amplitude is slow. One of the first implementation of this approximation was due to Caldeira and 
Leggett jl7f| for the problem of interacting quantum mechanical harmonic oscillators. Since then, generalizations of 
this kind of dynamics for the case of scalar fields have been implemented by several authors [|| -|l2|. In particular, 
in an earlier work by us Jl^j , a specific consistent solution regime for adiabatic dissipative dynamics was identified, 
in which the motion of the scalar field amplitude is overdamped. Our treatment in that work was restricted to the 
high-temperature regime. The purpose of this paper is to extend those results to zero temperature and to find valid 
parameter regions for the adiabatic dissipative regime. The key step is the inclusion of quasi-particles effects through 
the appropriate use of full two-point Green's functions at zero temperature. This plays a central role in the dynamics, 
as will be reviewed later in this paper. 

A primary, but not singular, motivation for understanding the overdamped regime is to realize the warm inflation 
scenarios from first principles. The extension of the overdamped solution to zero temperature has important conse- 
quences for this goal. A careful examination of the earlier work |l2| as well as (l^Jl^l reveal that the high-temperature 
overdamped solutions came very close to a full realization of warm inflation from simple interacting models. Although 
ultimately what we found for that regime was that full warm inflation solutions only could be constructed for certain 
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complicated models [|L8|-pO| , a foremost impediment to simpler solutions was the limitations of the high temperature 
approximation. The main accomplishment of this paper is to establish that the extension of the solution to T = 
can be adequately posed. Although only flat nonexpanding spacetime is treated in this paper, these results provide a 
necessary step towards the case of expanding spacetime. Furthermore, the results in this paper already are suggestive 
that once extension to expanding spacetime is done, warm inflation solution probably will be viable in most of the 
simple quantum field theory models. 

The calculation in this paper is incomplete in that dissipation at zero temperature necessarily implies a nonequilib- 
rium state, which immediately should be driven up to some excited statistical state. As such, the results presented here 
serve only as an indication for a nontrivial dissipative effect, which requires a more general nonequilibrium treatment 
for determining the precise nature of the excited statistical state. In this respect, the results in this paper are as yet 
insufficient for constructing completely consistent quantum field theory warm inflation solutions. Nevertheless, these 
calculations establish an important point, that interacting scalar field systems have an intrinsic and, in appropriate 
regimes, robust affinity to dissipate their energy. This fact in turn could place considerable doubt on a commonly 
followed, though unproven, assumption of inflationary cosmology, that dissipative effects of the scalar inflaton field can 
be ignored, thereby leading to a supercooled regime of inflation with vanishing, or negligible production of radiation. 

Despite the unconventional implications our results suggest for inflationary cosmology, the argument for supercooling 
in scalar field driven inflation deserves a reexamination for the following reason. The basic statement of supercooled 
inflation is that a single degree of freedom, the zero mode of the inflaton, maintains all the energy of the universe 
during the entire duration of inflation in the form of potential, or equivalently vacuum, energy. However, the curious 
point is, even if this single degree were to allow a minuscule fraction of the energy to be released, say one part in 10 20 , it 
still would constitute a significant radiation energy density component in the universe. For example, for inflation with 
vacuum energy at the GUT scale ~ 10 15_16 GeV, leaking one part in I0 20 of this energy into radiation corresponds 
to a temperature of I0 11 GeV, which is nonnegligible. In fact, the most relevant lower bound that cosmology places 
on the temperature after inflation comes from the success of hot Big-Bang nucleosynthesis, which thus requires the 
universe to be within the radiation dominated regime by T > lGeV. This limit can be met in the above example 
by dissipating as little as one part in I0 60 of the vacuum energy into radiation. Thus, from the perspective of both 
interacting field theory and basic notions of equipartition, it appears to be a highly tuned requirement of supercooled 
inflation to prohibit the inflaton from even such tiny amounts of dissipation. On the flip side, the warm inflation 
picture demonstrates pl| that by relaxing this requirement, the most unnatural and technically intractable aspect 
of supercooled inflation, reheating, becomes unnecessary. As such, despite the technical complications in computing 
dissipative effects, as exemplified by this and the earlier papers, on general grounds their presence should be a natural 
expectation. 

The paper is organized as follows. Our real-time dissipative formalism is reviewed in Sec. 0. Explicit expressions 
are given here for the real-time, fully dressed Bose and Fermi two-point Green's functions, which are valid in the 
entire temperature range, including T = 0. This formalism then is applied in Sects. || and [nj to obtain the effective 
equation of motion of the system, here a scalar field that has a classical amplitude out of equilibrium at T = 0, which 
is coupled to fields that act as a reservoir bath for dissipative energy exchange. Two types of models are treated 
in Sections JO] and [II, which we have denoted respectively as direct decay models, where the system field directly 
decays to light particles that are part of the environment, and indirect decay models, where this process is mediated 
by intermediate fields. In Sec. we apply the effective equations of motion to estimate the magnitude of radiation 
energy production. Finally in Sec. [v|, we give concluding remarks and comment on the implication of our results for 
the case of an expanding spacetime, which is relevant for inflationary cosmology. An Appendix is also included in 
which the computations are given of the zero temperature decay widths that are used in the paper. 



II. DIRECT DECAY MODELS 



In this section, models will be studied where the system, here a scalar field $, is coupled to bath fields of lighter 
mass than itself, thus allowing <!> to directly decay into the bath fields. In the first model the scalar field interacts 
with a set of fermion fields with Lagrangian density 

1 m 2 \ N * 

£[$, fo, fa] = -{d^f - - -$ 4 + Tpk [i - m^ k - h k $) fa . (2.1) 

fc=i 

The $ field is decomposed into a classical background component ip and a quantum fluctuation part as $ = <p + (f>, 
with 

($> = ¥>, (</>) = 0, 
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and where the background component is assumed to be homogeneous tp = (p(t). Using the tadpole methodf], where 
we impose (</>) = at all orders in perturbation theory, it leads to the condition that the sum of all tadpole terms 
vanish. The effective equation of motion (EOM) for ip(t) then becomes 



A A A N * 



(2.2) 



k=l 



where (4> 2 ), (</> 3 ) and (ipk^Pk) are given in terms of the coincidence limit of the (causal) two-point Green's functions 
G~^ + (x, x') and Si (x : x'). These Green's functions are obtained from the (1, l)-component of the real time matrix 
of full propagators which satisfy the appropriate Schwinger-Dyson equations (see, e.g., pX|Jl2[| for further details) 



□ + mi 



A 



G l f ) (x,x')+ / d zT, ( j ) (x, z)G < f,(z, x') = iS(x, x') 



and 



[i p- m^ k - h k ip] S^ k (x,x') + J d 4 z^ k (x, z)S^ k (z,x') = i5{x,x') . 
The momentum-space Fourier transform of G^,(x,x') (for the scalar field) can be expressed in the form 

G$(x, x) = i 



d 3 Q ai ^-^)fG+ + (q,t-t') G+-(q,t-t') 



(2tt) 3 



GT+(q,t-f) Gr( q ,t-t') 



(2.3) 



(2.4) 



(2.5) 



where 



G++(q, t-t') = G> (q, t - t')9(t - t') + G< (q, t - t')6(t' - t) 
G^-(q, t - t') = G> (q, t - t')0(t' -t)+G< (q, t - t')0(t - t'), 
G+-(d,t-t')=G<(d,t-t'), 
G- + (q,t-t>)=G>(q,t-t>). 



(2.6) 



In these expressions the fully dressed (field independent) two-point functions, at finite temperature T — 1//3, under 
the approximation that the spectral function for the scalar field has the standard Breit-Wigner form, are given by^| 

G>(q,< - t') = -L {[1 + n+(tJ4 - iT^e-^-^^-t") + n*(w* + iT^*^^ } 6{t - t') 

+ J- {[1 + n+fa + iT^}e-^+ ir ^ (t - t "> + - iT^y^-^^-^Xeit' - t) , 

G<(c l ,t-t')=G>(c l ,t'-t), (2.7) 

where is the Bose distribution function, uj^ = lj^{c\) is the particle's dispersion relation and is the </> decay 
width, defined as usual in terms of the field self-energy by 



lmE^q, w^) 
2clU 



(2.8) 



For the fcrmion fields we have instead that 








= s>( q ,t- 


t')6{t-t') - 


H5<(q,*- 


£')0(i' 


-*), 




*') 


= S£(cbt- 


t')9(t' - 1) - 


hS<(q,t- 


t')0(t- 


-0, 




*') 


= S$(q,t- 


0. 








s^ + (a,t- 


f) 


= ^(q,t- 


0. 









(2.9) 



1 For earlier references on the tadpole method applied to determining the equation of motion of a scalar field in the real time 
formalism, see for instance Ref. pi| . 
2 We thank Ian Lawrie for pointing out the correct form of these equations. 
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In this case S^' < are given by^| 

x [1 - n^luj^ - 

_ i(Uip+iT^){t—t') 



-i(uj^,—ir^)(t—t') 



1 



X [1 - n^{uj^ + iTi/,)] 



n^,[y>i> + «r.0) > 9(t - t') 



i(u$—iT$){t—t') 



-(City - ir^)7Q- ^ + m 4',r - « 



and 



2lo,i 



(2.10) 



x [1 - notify, + iT^,)] 



-i(a)w,— ir^,)(t— {') 



1 



,i{u)^,— t') 



x [1 - n^,(w^ - ?r^)] 



-(u;^ - ir^)7o— j4 + m i/>,r - * 



n^(u^ - iTy) > 0(t - t') 



(fcty + «r^,)7o- ^ + "ty,r + i 



n^(uj^ + iT^) > 0(t' - t) 



where — \ j(eP u * + 1), lu^ = y^q 2 + m 2 r , with m^ :J , = to^ + ReE^, and 



r^,(g) = ImE^o;,/,, q) . 



(2.11) 



(2.12) 



In this paper we will be interested in studying dissipation only of the background scalar field ip at T = 0. Therefore 
we will restrict ourselves just to the zero temperature expressions of the above equations. 

It will be assumed the couplings A, hk <C 1, so that perturbation theory can be co nsist ently formulated with 
subleading terms neglected. Then by perturbatively expanding the field averages in Eq. (|2.2| ), up to two-loop order 
in the scalar loops (higher order loop terms involving fermions also are higher order in the perturbation expansion, 
as we will discuss later), the effective EOM for ip(t) is 



( ^ + ^(0 2 ) O Jy(i) + ^ 3 (t) 



A 



+ A„-(/) / dt>^(t>) J -01m 



\2 ft 

-y / dtW) 



d 3 qi d 3 q 2 
(2tt) 3 (2tt) 3 



Irn 



G++(qi, t - t')G+ + (q_ 2 ,t - t')G++(qi +q 2 ,t-t') 



N,: 



fc=l 



/ dt' h k <p(t') I 7^x3 Im 



s+ + ( q ,t-t')s ++/3a (q,f'-*) 



(2.13) 



with (. . .)o meaning the vacuum expectation value at zero background field ip = 0. Using the explicit expressions for 
the Green's functions in Eq. (2.13) and at T = 0, this EOM becomes 



3 Analogous expressions for the general nonequilibrium case are given in 
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<p(t) 

3 



dt'<p(t') 



(2tt) 3 

d 3 qi d 3 q 2 sin[(^(qi) + ^(q 2 ) + c^(qi + q 2 ))|t - f'|] 



8w2 



xe 



(2tt) 3 (2tt) 3 

r*(qi)+r^(q2)+r^( qi +q 2 )]|t-i' 



-2 £ ft 

fc=l 



eft' <p(0 



8w (qi)u;0(q2)w <) !,(qi + q 2 ) 
d 3 g e - ar *»WI*-*'l 



(2tt) s 



{2I* fc w^ cos(2^Jt-t'|) 



[2mJ, r (wJ t - "4j + (c^ + m^ r - r^J] sia(2oty h |t - t'\)} = 



(2.14) 



The temporally nonlocal terms of the types appearing in the above equation have been shown in |2|,|11|- |14| , |24| , [25| to 
lead to dissipative dynamics in the EOMs. This can be made more explicit by an appropriate integration by parts in 
the time integrals in Eq. ( |2.14[ ) (see ||j2§]) to obtain the result (at T = 0) 

+ ml <p(t) + ~<p 3 (t) 

+\Mt) f dt'^t'Mt') [ /\ Kcos(2^|t-f|)+r,sin(2^|t-^|)] e . 2r4g)|t _, | 



(27T) 



8c 



+ Y / dt ' J ^j^|{K(qi) +^(q 2 ) +w (qi +q 2 )] 
x cos[(w0(qi) +w (q 2 ) +w (qi + q 2 ))|i-i'|] 

+ [r ( qi ) + r (q 2 ) + T (qi + q 2 )] sin[(w^(qx) + w (q 2 ) + w^(qi + q 2 ))|£ - 
e -[r*(qi)+r^(q2)+r^( qi +q 2 )]|t-t'| 



8cj0(qi)w0(q 2 )w0(qi + q 2 ) 



1 



[r (qi) + r (q 2 ) + r^(qi + q 2 )] 2 + [t^(qi) + ^(qfc) + ^(qi + q 2 )] 2 



fe=i 



dt' <p(t') 



,-2Tj, (q)\t-t'\ 



(27T) 



X {[r^^V-fc + 2t 4 fc m V>,r - 2 ™lr^ fc + ^fc^Vfc^* + m l,r)] COs(2w. Ji - t'|) 

+ [1^(2*4 - T% h )(u% h - m% >r ) - r^J sin(2^Jt - = , 



(2.15) 



where m$ and A are the effective mass and coupling constant given respectively by 



mA, + A 



1 



(2tt) 3 4w 
A 2 f d 3 qi d 3 q 2 



k=l 



d 3 q 

(2tt) 3 



o 



<4* 



(27r) 3 (2tt) 3 [ 8w </> (qi)w0(q 2 )u; <) !,(qi + q 2 )[w« ) i(qi) + w^(q 2 ) + u^(qi + q 2 )] 



w3 



(2.16) 



and 



A = A - A 



d 3 q 
(2tt) 3 



16^ 



■o(3 



(2.17) 



Both and A naively appear divergent in Eqs. ( [2.16] ) and (2.17) respectively due to the perturbative correction 
terms. However, they are rendered finite by the usual introduction of counterterms to renormalizc the mass and 
coupling constant in the original Lagrangian flljj ]. In this procedure, the bare (infinite) mass and scalar self-coupling 
in the Lagrangian are m — m^ r + Sm^ and A = A r + 5X, where the counterterms Srrij, and SX cancel the divergent 
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contributions in Eqs. (2.16) and (2.17) in the usual way. In what follows, we then can just interpret the masses and 
couplings appearing in our equations as the renormalized ones and we will omit any additional subscript and overbars 
for simplicity, so hereafter fh — > m and A — > A. 



A. The EOM in the adiabatic-Markovian approximation 



The effective equation of motion Eq. ( 2.11 ) is the main result of this section. As should appear evident, this 
equation is difficult to solve either analytically or numerically, since memory of the past history of ip is required 
at each stage of evolution. Falling short of a detailed analysis of this equation in this paper, we would like some 
r ough estimates of the dissipative effects described by the equation. In particular, as a first approximation to Eq. 



(2.15), we will explore the derivative expansion at leading order for the temporally nonlocal terms. This amounts to 



substituting if — ► t in the arguments of the fields entering in the time integrals. This is equivalent to a Markovian 



approximation for the dissipative kernels in Eq. (2.15). A solution regime for tp where such an approximation might 
be valid is the adiabatic regime, where the motion of ip is slow. In this and the next subsections, we will establish a 
self-consistent solution regime composed of both the adiabatic and Markovian approximations. To proceed, first the 
Markovian approximation will be implemented and then self-consistent solution regimes will be identified. Thus upon 
implementing the Markovian approximation, we then can easily perform the if integrals and obtain 

m + ml <p® + ^(t) + AV(W(t) / (0 8 ^4 + ^ )2 

A 2 f d 3 qi d 3 q 2 [r (qi) + I^(q 2 ) + r (q x + q 2 )] [w (qi) + w (q 2 ) + ^(q x + q 2 )] 



3 J (2tt) 3 (2tt) 3 4^(qi)cj0(q 2 )w0(qi + q 2 ) 

1 

X 2 

{[T^(qi) +I>(q 2 ) +I>(qi +q 2 )] 2 + [c^(qi) + w (q 2 ) + u; 4> (q 1 +q 2 )] 2 | 

+ £^"/p? \^k)- r *°°- (2 - 18) 

An additional simplification can be achieved by demanding that the (effective) field masses satis fy the ineq uality 
77i0 > 2m^ k . In this case, at zero temperature, the only contribution to the decay widths in Eqs. (2.8) and ( [2.12 ) come 



from the imaginary part of the one loop contribution to the self-energy for This is given by the internal fcrmionic 
propagators diagram, which then represents the decay rate for the kinematically allowed process <& — > ipi- + We 
therefore have that T^ k = and r is given by (see the Appendix) 

N * U2 ( Am 2 \ ^ 



We must point out that even though in this case the contribution to the dissipative term in Eq. ( 2.18| ) coming from 
the fermionic loop vanishes, i.e. T^ k = 0, this is a consequence of assuming the Mark ovian approximation for the 
dissipative kernels. In the general case, the last term in the non- local EOM, Eq. ( §If ), with I> fc = 0, still can be 
interpreted as a dissipative contribution |^|l3|,[l4| to the EOM. Therefore, h ere t he Markovian approximation somewhat 



under-estimates the whole dissipative nature of the non-local kernels in ( 2.15 ). We expect this not to invalidate our 
main objective here, which is to determine the viability of adiabatic dissipative behavior in the background field 
evolution and the intrinsic dissipative nature of the field dynamics. 

As an aside, it is inte restin g to discuss briefly the interpretation of dissipation coming from the second non-local 



dissipative term in Eq. (2.15) when = 0. Though for this case, particle decay does not contribute to dissipation, 
there still is dissipation due to "off-shell" excitation of virtual states which leads to decoherence and power law decay 
of the background field a mplit ude 0], in the absence of the scalar field self-interactions. This contrasts with the first 



dissipative term in Eq. ( [2.15 ), in which dissipation truly is coming from the real scattering by the quasi-particles 



in the medium. Furthermore, dissipative kernels of the sort treated here typically have some type of long time tail, 
which retains memory from the past. However, as will be discussed later, a non-vanishing decay width, as present in 
this case, helps to suppress the long time tail of the kernel. 

Studies of terms similar to the second nonlocal term in Eq. ( [2.15| ) have been done in [^6] . These works studied 
the case N$ — 1 and T$ — 0, and examined the linearized form of the effective EOM once the fermion fields were 



G 



integrated out. An important issue raised in p7j , p8[ was in regards the s ingu larities in the EOM at the initial time, 
which in our case refers to the behavior of the dissipative kernel in Eq. ( 2.15 ) arising from the fermionic loop, when 
computed at t = t'. Observe that this term appears divergent in the ultraviolet region, but this apparent problem 
easily can be solved in our approach. To see that our final equation in this case is ide ntical to the final result obtained 
by the authors in [ p7|pg| ], take for example the second nonlocal term in Eq. ( [2.15 ), due to the fermionic loop, and 
integrate the time integral twice by parts to give 



■4> 

E 

k=l 



dt' ip(t') 



d 3 c 



-2IW (q)\t-t'\ 



(2tt) 3 



x { \F% k ^ k + 2uj l k m % k ~ 2m % k ^ k + v % k ^ k {^% k + m |j] cos(2w^Jt - t'|) 
+ [r^(2m^ - I^J(w^ -m\ h )- I* J sin(2^Jt - t'\)} 

d 3 q ( l 




+3-^ sin(2w^ k \t - *' 



(2.20) 



The first term on the rhs of ( 2.20| ) is just the second dissipative term appearing in Eq. (2.18). By taking Y. l p k — this 
term vanishes. The second term is a wave-function renormalization term, which can be absorbed in a wave function 
counterter m (SZ ) in the Lagrangian, by rewriting the kinetic term as (<9 M $) 2 — > (1 + 5 Z)(d u & ) 2 . The last term on the 
rhs of Eq. J2.20|) , for T^ k = and N$ — 1, reproduces exactly the result obtained in |p7| , p8| for the EOM integrated 
over the fermion fields. Thus, we sec that there is no ambiguities or diver genc es here associated with the kernels 
evaluated at equal times. Also, the non-local two-loop scalar term in Eq. ( 2.15 ) can be examined in a similar way 
and from such an analysis its contribution to the wave- function renorma lizati on can be extracted. 

Returning from this digression to the EOM Eq. ( 2. IS ), and using Eq. (fyH) for T and I> = 0, the EOM becomes 



0(f) + m\ <p(t) + -^ 3 (f) + V(<p)<p(t) = • 



(2.21) 



Here T)(ip) is the dissipative coefficients, which after doing the momentum integral in the "one-loop" dissipation term 
and using the symmetry of the "two-loop" dissipation term under change of momentum integration variables, becomes 



X 2 a 2 




4> T 

d 3 qi d 3 q 2 f 



2m' 



(2ir) 3 (2ir) 3 1 W0 (q 1 ) 2 W 0(q 2 )w0(qi + q 2 ) [c^(qi) + ^(q 2 ) + ^(qi + q 2 ) 



(2.22) 



with 




(2.23) 



Finally, we also can easily work out an equivalent model of $ coupled to bath fields that rather than fermionic fields 



are scalar fields, Xji 3 = 1 • • ■ A^ x , with a trilinear coupling as Xw=i if^xf • Once again, by choosing m^ > 2m Xj we 
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have T x . = but a nonvanishing $ de cay w idth, T^, which also has been evaluated in the Appendix. Going through 



the same steps as used to obtain Eq. ( 2.21 ), we obtain an analogous expression with a similar dissipative coefficient 
r](tp), except with — > a,p. x , where 




Am 



(2.24) 



B. Examination of the dissipative kernels 



Up to this point 
toEq 



Eq 
Eq. 
general forr 



( 2.15 
2.15 



the naive implementation of the Markovian approximation has been examined and has led from 
The applicability of this approximation is now considered. Returning to the no n-loca l EOM 

in the 



(2.21 



we can express for example the term generating the one-loop dissipation contribution in Eq. (2.21 



dt' tp(t')tp(t')K(t,t') , 



where the dissipative kernel is given by 



K(t, 



t') = X 2 [ ^ cos ( 2uj <p\t ~ + r sin(2i^,|t - t'\)} ^ 2 t s ( q )\t-t' \ 



(27T) 



(2.25) 



In Fig. 1 K(t,t' — 0)/A 2 is plotted for four cases, T^{Q)/m^ = 0.1,0.5,1.0,5.0. There is a degeneracy amongst the 
parameters Ny, h and m.0 from which to choose for a given value of 1^(0). The plots of the kernel are given over 
four different time interval, in order to see the different aspects of its behavior. 

As can be seen, the kernel has a pronounced, narrow peak around t = t' followed by a power-law de cayin g oscillatory 
behavior, with greater decaying as the nu mber of fields coupled to <I> increases. Since from Eq. (§j||) r ~ l/|q| 
for |q| 3> it can be seen from Eq. ( p. 25 ) that once |q| > 2m^T(0)\t — t'\, damping from the exponential 
term relinquishes. Thus for \t — t'\ > l/m$ the behavior of the kernel is K(t,t') ~ — X 2 Ci(4 : 'm r j > T r j 1 (0)\t — t'\ 2 ) ~ 
-A 2 sin(4m^r^(0)|t - t'\ 2 )/ '(4m^r^(0)|t - t'\ 2 ). Note this long time behavior of the T = contribution of the kernel 
differs from its high temperature component which was studied in [|ll 25 1. At high temperature, the decay widths at 
large |q| behave as ~ T 2 /|q|, but there are also factors of the number density which become Boltzmann suppressed. 
As such, the high temperature limit of the kernel becomes highly exponentially suppressed at la rge t imes. 

Power-law decay of the kernel at T = implies that memory of the scalar field is retained in Eq. ( 2.1£ ) in determining 
its future evolution. As such, the derivative expansion clearly is not generally valid. However, for sufficiently slow 
motion of the scalar field, the derivative expansion still may be valid for some duration of time. In particular, suppose 
tp/ tp sa 7 in a given solution, where 7 is the approximate magnitude of this ratio over some interval of time. Then, 
self-consistency of this solution with respect to the derivative expansion holds for a time interval \t — t'\ < 1/7. Thus, 
as mentioned at the beginning of this subsection, the slower is the motion of tp, the longer the above approximation 
is valid. Such slowly varying solutions are useful to investigate due to their simplicity. They also may have practical 
use for example for warm in flatio n, where one seeks solutions where the motion of tp is slow. 

Although the kernel Eq. ( 2.25Q does retain past memory, it is worth noting that on general grounds for reasonable 
tp(t) solution regimes, the memory retention only is up to some finite time in the past. The observation here is that 
the oscillation rate of the kernel increases with time due to the quadratic dependence on time K(t, t') ~ sin(a|f — 1'| 2 ), 
where a is a constant. As such, for reasonable motions of tp, beyond some time interval into the past \t — t'\ > Aio, 
the characteristic oscillation frequency of the kernel will exceed that of tp. Thus at all times past this point, the 
contributions from tp(t) primarily cancel. As a practical point, despite this property of the kernel to filter through 
increasingly high frequency components as the time passes, there are limitations to the types of tp(t) motions which 
it can describe. In particular, this feature of the kernel most efficiently is able to cut-off the long time memory if 



* Although here we analyze the one-loop dissipativ e ker nel in detail, similar conclusions also can be shown to apply for the 
more complicate two-loop dissipative kernel in Eq. (2.15). 
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the field configuration only has slow frequency components. As such, ideally the applicability of this kernel is near 
equilibrium conditions. Note, this also is the regime assumed in previous works [fll| , |l2[ for the finite temperature case. 

An additional consistency condition related to the adiabatic approximation at finite T, considered in jl2],[l8]], is 
that the microscopic dynamics should be much faster than all macroscopic motions. For the finite-T case, this was 
a necessary requirement since it guaranteed the system thermalized fast enough to adjust to any changes in the 
macroscopic state. The need for rapid thermalization was necessary for self-consistency of the solutions. In contrast, 
the T — dynamics treated in this paper does not require a specific statistical state in which the scalar field evolution 
occurs. Thus, consistency requirements, if any, with respect to the rates of microscopic versus macroscopic dynamics 
are less well defined for the T = case. To properly address this question, a complete nonequilibrium analysis 
is necessary, which is beyond the scope of this paper. However, for the time being, to be conservative, a similar 
consistency condition to the finite-T case adopted in p^ j will be imposed, which requires the rate of microscopic 
physics to be faster than all macroscopic motions. Since the only scales characterizing the microscopic physics are 
the decay rates, this requirement implies 



> r _1 , (2.26) 



which is analogous to the condition in ]12{ , except above the zero temperature decay rates are used. In Sec. IV some 
estimates will be given of dissipative dynamics based on the combined adiabatic-Markovian approximation of this and 
the previous subsections. 



III. INDIRECT DECAY MODELS 



This section will consider models in which the particles ultimately created from dissipation of the scalar field are 
coupled indirectly to the scalar field through an intermediate field. Such a case has much more variety in the types of 
decay sequences, as compared to the direct decay models of Sec. EL The basic model to be examined consists of the 



system, a scalar field $, along with scalar fields %ji 3 = 1 
density is given by 



. N x and fcrmion fields ipk, k = 1 . . . N^. The Lagrangian 



1 in \ 



N, L 



k=l 



2 m X 3 - ,1 



i ft - m^, k - h k ,<j,® - ^2 h kj, x X] 



2 3 4! 3 2 3 



ipk , 



(3.1) 



where all coupling constants are positive: A, fj,g 2 , hk,<j>, ^&i,x > 0. 

As before, we are interested in obtaining the EOM for a scalar field configuration cp = ($). For this, the fields Xj 
and tpk are regarded as part of the environment. Once again the scalar field $ is decomposed into its expectation 
value and fluctuation, $ = </? + </>, where ($) = ip. The EOM for cp then is obtained from the tadpole method by 
imposing that ((f)) =0, which leads to the condition that the sum of all tadpole terms vanish. Restricting again our 
analysis of the EOM to a homogeneous field ip = ip(t), we obtain the effective EOM for tp 



(p{t) + m%<p(t) + ^(t) + ^ip{t)(<jP) 



X 



J'=l 



a) [<p(t){x. 



A', , 



■/]hk,<t>{ipkj>k) = , 



(3.2) 



fc=i 



where the field averages above can be expressed as usual in terms of the coincidence limit of the (causal) two-point 
Green's functions G^ + (x, x'), G^ + (x, %') and S^ + (x, x') for the $, Xj an d ^fc fields, respectively. Using the expressions 
of the previous Section for G ++ a nd S ++ for the scalar and fermionic propagators respectively, and working out the 
expression analogous to Eq. ( 2.15| ), we obtain the EOM 



<p{t)+ml (p(t) + ^<p 3 (t) 
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+ Y / / ^i^^{[^( c n) + ^(q2) + ^(qi + q2)] 

x cos[(o^(qi) + o>0(q 2 ) + w^(qi + q 2 ))K - i'| ] 

+ [r (qi) + r (q 2 ) + I>(qi + q 2 )] sin[(^(qi) + w (q 2 ) + u^(qi + q 2 ))|i - t'\ }} 

e -[r*(qi)+r^(q 2 )+r^( qi +q 2 )]|t-t' 



8w0(qi)w0(q 2 )w0(qi +q 2 ) 

1 



x 



[r^(qi) + iyq 2 ) + T (qi + q 2 )] 2 + [o^(qi) + w (q 2 ) + w^(qi + q 2 )] 2 



+ X^ 4 / tft 7 ¥>(0 y ^I^^ {K,(qi) +^(Q2) +^ (qi +q 2 )] 
x cos[(w Xj (qi) +w (q 2 ) +w (qi + q 2 ))|i-i'|] 

+ [T Xj ( qi ) + r^(q 2 ) + r (qj + q 2 )] sm[(uj xj ( qi ) + ^ (q 2 ) + w (qi + q 2 ))|i - t'\ ]} 

e -[r Xj (qi)+r^(q 2 )+r^,( qi +q 2 )]|t-t'| 
v 

8w Xj (qi)cj (q 2 V (q 1 + q 2 ) 

1 

2 2" 

[r x ,(qi) + T (q 2 ) +r^(qx + q 2 )] + [uj Xj (q x ) + w (q 2 ) + w (qi + q 2 )] 



■XX, / dtw) / 



d 3 q e -2r^ k ( q )\t-t'\ 



x { [ r V' fc ^V' fc + 2uj l k m l k ~ 2n A )k u)^ Jk + T^ h w^ k {ujl k + m|J] cos(2w v , fc |f - 

+ [r^ fc (2m2 fe - r 2 _ m 2 J - I*J sin(2c^Ji - t'\)} = . (3.3) 

The dissipative dynamics in these indirect decay models will differ based on the relation of the masses amongst the 
<&, Xj an( i i'i fields. Consider first the case where the (effective) masses satisfy the relation 

> m Xj > 2m^ k , (3.4) 

which we will refer to more specifically as the indirect cascade decay regime. The first attribute of this regime to 
note is at zero temperature and T Xj are non-vanishing and T.^ — 0. There are two kinematically allowed on-shell 
processes, $ — > ipk + ipk 1 with decay width as given in Eq. ( 2.19 ), and Xj ~~ * ipk + V'fcj with decay width 



N,p h 2 ( Am 2 * 



r ( q ) = yJ^ml. 1 . (3.5) 

Implementing once again the adiabatic-Markovian approximation of Subsec. IIA, we then obtain the same expression 
for the EOM as Eq. (2.21), except rj((p) now is given by 

\ 2 2 

77(99) = if (t) 



X 2 »7..,. f <f\; ,{\- 2 J 1 +o(^ 

(2tt) 3 (2tt) 3 ) ^( qi )2 W0 (q 2 )^( qi + q 2 ) [^(qi) + w (q 2 ) + w,(qi + q 2 )] 3 I u\ 
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2 

4 a x> 



3'=1 



32tt 



r/v 



x 3 - 



Ax sl a x,v> f d3 qi d3 q2 



E 

-o 



(2tt) 3 (2tt) 2 



'r 2 r 2 



l x,V> 



2.,/m 



2m 2 



; (qi) 2 w (q 2 )w0(qi +q 2 ) [w Xj (qi) + w (q 2 ) + w«^(qi +q 2 )]' 



(3.6) 



with a^j given by Eq. (2.23) and 



x,V" 



N 4- h 2 

Ekj,x 2 
— — mz . 

k=l 



(3.7) 



It is interesting to note that the last two terms in Eq. (3.6) come solely from the decay channels of Xj m to the 
fermion fields ipu, which then backreact on the system field tp in terms of a damping force. The same effect also would 
appear for the case where $ was the lightest field, 



i Xj > 2m^ k > 7710. 



(3.8) 



This regime will be referred to simply as the indirect decay regime outside the cascade region Eq. (3.4). For this case 
= 0, or a^tf = in Eq. ( |3.6[ ) . This type of dissipation, in which the system field is lighter than the decay products 
was first noted by Calzetta and Hu in |14 1 . They have shown how a heavy field influences the dynamics of a light field 
in the form of dissipation and fluctuations of the light field, even when no aspect of the light field dynamics is above 
the mass threshold of the heavy field. This same behavior also can be inferred from our results for the above regime. 

Let us finally make a few c omm ents abou t higher order loop terms. Observe that the results obtained for the 
dissipation coefficients in Eqs. ( 2.22 ) and ( |3.6D are the leading order ones at zero temperature, 0{)?h\) and 0{gjh%j X ) 
respectively. We have neglected higher order loop contributions to the EOM, since they all can be shown also to be of 
higher order in the coupling constants. For instance, a two-loop contribution made of a fermion loop with a vertical 
scalar propagator can easily be seen to give a contribution to the dissipation coefficient r](ip) of order 0(h^T^,) ~ 
for the first case analyzed in Sec. [n], and 0{h\ ), for the case studied in this section. Also, higher order scalar 

loop terms are subleading in the coupling constants, as compared to the results given by Eqs. ( 2.22 ) and ([T(]). 
For example, at finite temperature, scalar ladder diagrams Jl^,|2^] are known to be of the same order as the two- 
vertex one-loop term in the EOM, due to on-shell divergences of these diagrams and the way the field decay widths 
regularize them. However at zero temperature these divergences are not present and ladder diagrams are at most of 
order C(A 4 r 2 ) - C(A 4 ^), for the model in Sec. || and 0(g*h 4 kj x ), for the model in this section. 



IV. APPLICATION 



In this Section we will apply the effective equation of motion Eq. (2.21), with rj(ip) as given by Eqs. (2.22) and (|3.6|), 
to make estimates of entropy production from conversion of the scalar field potential energy into radiation. Recall the 

fplicability of Eq. (2.21) is limited, since it requires the validity of the adiabatic-Markovian approximation of Sect. 
For this reason and since as mentioned in the Introduction, the zero temperature dissipation found here is indicative 
of a nonequilibrium dynamics that drives the system to finite temperature, we will not delve into detailed applications 
of this effective equation of motion. The full dynamics of this problem must be understood before detailed application 
is worthwhile. Nevertheless, we will make some naive estimates of radiation production from our equations, just to 
get a feeling for the magnitudes of the effect. Note that for application to warm inflation, the fact dissipation occurs 
at T = implies the dynamics automatically will bootstrap the universe to finite temperature, independent of initial 
conditions. 



In general for a damped equation of motion of the form (2.21) the regime of overdamped (underdamped) motion is 



(4.1) 



For the direct decay models in Sec. [H], since the dissipative coefficient rj(ip) is suppressed by two powers of the <f> self- 
coupling parameter A, we find that underdamped motion generally is possible, unless there are very many environment 
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bath fields. For the indirect decay models of Sect. III. the cascade decay regime Eq. (3.4) has similar co nstra ints, 
thus solutions, as those of the direct decay models of Sect. ||. However, in the other interesting regime Eq. (^8), the 
dissipative coefficient and the parameters of the ip potential can be independently tuned, which means the overdamped 
regime can be identified for as few as one or two heat bath fields. Both the underdamped and overdamped regimes 
may have application to the chiral transition in heavy ion collision E3], whereas the overdamped regime also is of 
interest to the warm inflation scenario Below, both the direct and indirect de cay models will b e ex amined. 
For the direct decay model, the overdamped (underdamped) regime from Eqs. (2.22), ( 2.23| ) and (4.1) is given by 



A 2 AV 
2048tt 5 



> (<)1, 



(4.2) 



where the Yukawa coupling h is determined by requiring ~ hip < m(ip)/2 w y/Xtp/2. Thus, even for strong 

coupling A ~ 1, it requires > 10 fields for the overdamped regime. 

Considering first the underdamped regime, the adiabatic condition Eq. ( |2.26| ) requires 



ip 

— w m(ip) < Fj, 



h 2 N^ 
8tt 



i(cp) 



16tt 



■m(ip), 



(4.3) 



where for underdamped motion from Eq. (2.21) ip < m(ip)ip. The energy dissipated by the scalar field goes into 
radiation energy density p r , here composed of fermions and/or scalar bosons, at the rate 



Pr 



dt 



T)(<p)(p 



A 2 ^m 5 (y) 
2048tt 2 ' 



(4.4) 



where to obtain the last expression we estimate ip m m (tp)<p . For the overdamped regime from Eq. (2.21) 



(ip)(p/r)((p). Writing r)(ip) = Qm(ip), from Eqs. ( |2.22[ ) and ( f2~25| ) Q « A 2 A^/(2048tt 2 ) and for overdam ping E q 



(4.1) requires Q > 1. This requirement of overdamping automatically implies the adiabatic condition Eq. ( 2.26| ) is 
satisfied. Thus the radiation production in this case is 



m 3 ((p)ip 2 



(4.5) 



Overall, for the direct decay model, the underdamped regime is generic except if there are a very large number of 
bath fields N > 10 4 , in which case overdamped motion also becomes possible. In the underdamped regime, moderate 
radiation energy production occurs. In particular during a characteristic oscillation time ~ l/m(ip), the produced 
radiation has an associated temperature scale T ~ (pr/N^p) 1 ^ 4 w A 1 / 2 m((/?)/10 < m(ip). On the other hand, the 
overdamped regime, although requiring a large number of bath fields, can yield sizable radiation by increasing the 
amplitude ip in Eq. (4.5) 

Turning next to the indirect decay models of Sec. 



Ill, the cascade region Eq. (3.4) leads to similar solutions as 



given above for the direct de cay models, so will not be further elaborated. However, dissipative dynamics also appears 
to occur in the regime Eq. (3 



which does not have a direct interpretation in terms of particle decay at one of the 
two steps of the process, the ip — > \ transition. We believe further investigation is needed of this case in order to 
obtain a sensible interpretation of this process, and this is left for future work. Nevertheless, here it is interesting to 
estimate the size of the dissipative effects for this case. In particular for this case, the overda mpe d regime can easily 
be obtained as will be shown next. The overdamped regime requires the condition in Eq. (4.1) and the adiabatic 
condition requires 



m 2 (ip) 



(4.6) 



Since the scalar field sector has two free parameters, and A, it is always possible to tune m 2 (ip) to satisfy both 
the above requirements independent of T x , r](ip), and for any amplitude tp. 

In the overdamped regime, the kinetic energy of the scalar field is negligible. Thus, the loss in its potential energy 
translates into the energy released into radiation p r as 



p r (t) = r](ip)4> 2 



dV 

dp 



(mlip+-ip 3 )ip^V(ip)'- 



< V(p)T x , 



(4.7) 
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where the second line follows from Eqs. (2.21) and (3.1). Generally the microscopic scale is determined by T x . In 
this time interval, we find the radiation to increase to 



Pr {l/T x )^V{ V )^^<V{ V ), (4.8) 
v 1 x 



where to obtain the right most expression we used Eq. (1J3). Thus the energy dissipated into radiation is proportional 



to the potential energy contained in the scalar field. To consider some numbers, for example typical for inflation, 
suppose the potential energy is at the GUT scale V(ip) 1 ^ 4 ~ 10 15_16 GeV and m 2 {ip) / (r]T x ) w 10~ 4 . For this, it implies 
a radiation component is generated which, if expressed in terms of temperature, is at the scale T ~ 10 14 ~ 15 GeV, and 
this is nonnegligible. The main point to note is that considerable radiation production can occur from the indirect 
decay models in the regime Eq. (|3.8|). 



V. CONCLUSION 



Although the calculations in this paper were for non-expanding Minkowski spacetime, a brief reflection will be made 
here on the consequences of these results for inflationary cosmology. In this case, the background space is expanding. 
The effect of the interaction between the scalar field, in this case called the inflaton, and the background metric is 
known to yield a "SHip term in the inflaton effective equation of motion, where the Hubble parameter H = a/a and 
a{t) is the cosmic scale factor. As well known, this term does not arise from microscopic interactions with other fields, 
but rather from the macroscopic interaction with the background metric of gravity. Precisely this disparity in scales 
and the difference in origin of the interactions suggests that the effect of this term and the dissipative term computed 
in this paper will act independently on the inflaton to a good approximation, with perhaps some self-consistency 
requirements. At present, we are extending our calculation to expanding spacetime in order to examine this point. 
Should this work render true the expectations from Sec. IV for dissipation, it will be difficult to justify the supercooled 
inflation picture and rather it would appear the warm inflation picture is the natural one. 

For the moment, assuming the correctness of the above expectations, a lower bound on the temperature of the 
universe during inflation can be estimated from the direct decay model. In typical inflation models, m$ ~ H for 
supercooled inflation and perhaps a few orders of magnitude bigger than H for warm inflation. Furthermore, the 
conditions on density perturbations generally require the inflaton self-coupling parameter to be tiny A ~ 10~( 10 ~ 16 ). 
In a n exp anding background, the evolution of the radiation in presence of a dissipating scalar field source is given by 
Eq. (4.3) with the addition of the term —AHp r to the right hand side, which accounts for the redshift of the radiation 
due to background expansion. Therefore, in steady state p r — 0, p r w ijip 2 /{AH), or the associated temperature 
T ~ {pr/g*) 1 ^, where we will take the number of light particles g* ~ N. With these estimates, based on Eq. (|4.3j ) 
we find T > X 1 / 2 m 4> /10. If m - 10 10 GeV as a typical value, then this implies T ~ (10 4 - 10)GeV. 

In order not to affect the successful predictions of nucleosynthesis, the primary requirement is that the universe 
should be well within the hot Big-Bang radiation dominated regime by T « lOMeV. Slightly more conservative, 
though not necessarily mandatory, is to require that the QCD phase transition at T < lGeV occurs within the 
radiation dominated era. So a safe lower bound for inflation to end and the radiation dominated era to commence 
is T > lGeV. As such, the lower limits for radiation production during inflation given above still would be above 
this lower limit requirement set by cosmology. The results found here may be also useful in applications to the low 
temperature regimes of warm inflation identified in the phenomenological studies of [pl| . 

It should be clarified that the results found in this paper in no way require supersymmetry. However, these 
calculations easily could be applied in SUSY models. In such models, it is becoming appreciated that to avoid 
gravitino overproduction, for any type of inflation scenario the temperature of the universe after inflation can not be 
very high T < 10 10 GeV [ [32[ . From this perspective, the possibility found in this paper for low temperature warm 
inflation solutions in the direct decay models would be phenomenologically attractive. 

One cautionary remark is in order. Since ~ H, the lower limit on the temperature during inflation that is 
suggested above, is below the so called Gibbons-Hawking "temperature" Tqh = H/{2tt). However Tqh does not 
represent a temperature in the usual sense of a thermal bath of particles. Tqh acts like a temperature in the formal 
sense that for a non-interacting scalar field in de Sitter space, its euclideanized de Sitter invariant Green's function 
is periodic in imaginary time. The role that Tqh plays in this Green's function is formally the same as what actual 
temperature plays in the static thermal Green's function in Minkowski space. However, for the non-interacting scalar 
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field in the de Sitter case, there are no particles present in the sense that the field is in a vacuum state[|. In contrast, the 
radiation production we have computed for the interacting scalar field results in real particle production, irrespective 
whether its associated temperature is above or below Tqh ■ 

In summary, this paper has studied dissipative effects of interacting scalar fields at zero temperature. Similar 
treatments along this line are limited |^|,^|,^| and one of the novel features of this paper is the appropriate inclusion 
of quasi-particlcs effects through the fully dressed zero temperature two-point Green's functions. Another feature of 
our analysis which has been studied only to a limited extent in the literature is a detailed examination of the dissipative 



kernel in Subsec. [IB. The models examined in this paper were generic and in all cases dissipation was found. Since 
dissipative effects are seen for the zero temperature state, we conclude that radiation production from dissipation is 
invariably present for generic interacting scalar field systems, although the extent of radiation production can vary 
immensely. Minimally, it appears the mass of the scalar field times a suitable dimensionless coupling constant sets 
a lower limit to the associated temperature scale of the produced radiation. However, for the indirect decay models 



in the region Eq. (3.8), there appears a much more robust possibility for producing radiation. Although formally 
this is what is indicated by our calculations, as mentioned earlier, we feel further investigation is necessary of these 
indirect decay models in order to obtain a sensible interpretation of its dissipative process. In regards the potential 
implications of the results found in this paper to inflationary cosmology, we infer that under generic circumstances the 
scalar inflaton field will dissipate a nonnegligible amount of radiation during inflation. In particular, the lower bound 
suggested by the above estimates for the direct decay models already are sufficiently high to preclude a manda tory 



requirement for a reheating. Furthermore, the upper bound from the indirect decay models in the regime Eq. (3.8) 



could yield very high temperature warm inflation solutions, in the range discussed below Eq. (4.8). However, these 
only are expectations suggested by the calculations in this paper. Verification of these expectations requires a proper 
extension of these calculations to expanding spacetime, which we currently are examining. 
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APPENDIX A: 

In this appendix the decay widths for the processes 4> — > ip + (or x ~ * V* + w) arL d 4> ~~ * X + X are derived. Recall 
the basic expression for the decay of an initial particle of momentum p into two particles, 



r(p) = ^~ 

2w p 



X S(kf - to? )5{kl - m 2 ) J] (2m.,) ^ \M 



M 2 



fcrmion. 



(Al) 



where ui p — y^p 2 + M 2 , M is the mass of the scalar decaying field, and m stands for the mass of the decay products. 
Observe that the expression in the square brackets is Lorentz invariant, thus most conveniently it is evaluated in the 
rest frame of the initial particle p = 0. 

For the scalar to 2 fermion model [h is the Yukawa coupling) 

d = -h<j>(x)$(x)ip(x) , (A2) 

which implies 

M = -ihu sl (ki)u B2 (k 2 ), (A3) 



5 We thank Larry Ford for this clarification 
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so that 



J2\M\ 2 = k 



, ki ■ k% — m 2 



(A4) 



Substituting this into Eq. (Al) gives 



r(p) = 



h 2 M 2 ( Am 2 



2\ 3/2 



1 



87rw(p) V M 2 

For the scalar to two scalar model, with coupling constant g, we have instead that 



(A5) 



which implies 



£/ = -~2<l>{x)X 2 {x), 



M = -ig 2 



(A6) 



(A7) 



and so 



rr )= g4 (\ 4m2xV2 

W 16^(p) ^ M 2 



(A8) 
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Figure Caption 



Fig. 1. The kernel K(t,0)/X 2 in Eq. (^25|) for the cases r (O)/m 
5.0 (dot-dashed), plotted over four different time intervals. 
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